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Abstract 

Microphone directional array technology 
continues to be a critical part of the overall 
instrumentation suite for experimental 
aeroacoustics. Unfortunately, high sensor cost 
remains one of the limiting factors in the 
construction of very high-density arrays (i.e., 
arrays containing several hundred channels or 
more) which could be used to implement 
advanced beamforming algorithms. In an effort 
to reduce the implementation cost of such arrays, 
the authors have undertaken a systematic 
performance analysis of a prototype 
3 5 -microphone array populated with commercial 
electret condenser microphones. An ensemble of 
microphones coupling commercially available 
electret cartridges with passive signal 
conditioning circuitry was fabricated for use with 
the Langley Large Aperture Directional Array 
(LADA). A performance analysis consisting of 
three phases was then performed: 
(1) characterize the acoustic response of the 
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microphones via laboratory testing and 
calibration, (2) evaluate the beamforming 
capability of the electret-based LADA using a 
series of independently controlled point sources 
in an anechoic environment, and (3) demonstrate 
the utility of an electret-based directional array in 
a real-world application, in this case a cold flow 
jet operating at high subsonic velocities. The 
results of the investigation revealed a 
microphone frequency response suitable for 
directional array use over a range of 
250 Hz - 40 kHz, a successful beamforming 
evaluation using the electret-populated LADA to 
measure simple point sources at frequencies up 
to 20 kHz, and a successful demonstration using 
the array to measure noise generated by the cold 
flow jet. This paper presents an overview of the 
tests conducted along with sample data obtained 
from those tests. 

Nomenclature 

Dj free expansion jet nozzle 

diameter, ft 

e array steering vector, 

see eqn. (4) 

/ frequency, Hz 

G cross spectral matrix 

k acoustic wave number, ft' 1 

K number of steering locations 

M total number of array 

microphones 

N number of data blocks 


1 

American Institute of Aeronautics and Astronautics 



Pie) 

array output power for steering 
vector e , see eqn. (5) 

Pf 

RMS pressure of fundamental 
plus harmonic frequencies, Pa 

Ph 

RMS pressure of harmonic 
frequencies, Pa 

St 

Strouhal number 

r° r° 

distance from source to X c 
and m th microphone, ft 

r ^,n 

distance from steering location 
to x c and m th microphone, ft 

R k 

normalized array output of 
reference LADA at k th 
steering location 

THD 

total harmonic distortion, 
percent, see eqn. (1) 

T k 

normalized array output of 
ECM LADA at k th 
steering location 

Uj 

free expansion jet exit velocity, 
ft/sec 

W(k,x,x°) 

ideal array response, 
see eqn. (7) 

W, 

hamming window weighting 

X 

Cartesian coordinate, ft 


array phase center, ft 

X ik (f) 

k th FFT data block for i th 
microphone, see eqn. (3a) 

x m 

distance vector from steering 
location to m th microphone, ft 

^ RMS 

weighted RMS array response 
error 

CO 

angular frequency ( = 271 f), 
radians / second 

Introduction 


Aeroacoustic testing capability in both 
anechoic and hard-walled facilities has grown 
tremendously over the last several years, with 
successful collection of noise source location and 
directivity data in a variety of experiments. One 
of the reasons for the success of these 
measurements has been the development of both 
microphone directional array (DA) 
instrumentation and associated robust data 
processing algorithms applied in the analysis of 
the acquired data. 1 ' 7 While the gain in test 
productivity resulting from the use of DA 
systems has been significant, several challenges 
remain to be overcome before fully realizing the 


potential of this technology. For instance, the 
next generation of array processing algorithms 
being developed will potentially require the 
acquisition of significantly more acoustic 
information than is presently done, in some cases 
requiring acquisition of acoustic boundary 
conditions over the entire sphere of source 
directivity. DA systems involved in these types 
of tests will employ hundreds or thousands of 
individual microphones. Traditional array 
hardware, typically composed of externally- 
biased condenser microphone / preamplifier / 
power supply combinations for a number of 
acquisition channels, is in many instances too 
expensive to implement a high-density 
microphone system. Therefore, it is desirable to 
identify devices and hardware which can be used 
to implement a system exhibiting measurement 
performance comparable to the current 
generation of arrays, but at much lower cost. 
Several alternate microphone devices have been 
described in the literature which exhibit 
characteristics suitable for DA use and at 
economical cost. Among these are MEMS- 
based devices described by Sheplak 8 , Saini 9 , and 
Arnold 10 , and thin-film electret devices described 
byHsieh et.al. 11 

Over the last few years, the rapid pace of 
micro-fabrication technology development has 
yielded several promising new commercial 
microphones. Given the increasing interest at 
Langley Research Center and elsewhere in 
building higher density array systems, the 
authors initiated a systematic search for a 
suitable low-cost commercial off-the-shelf 
(COTS) microphone which could be used for 
this purpose. After examining several types of 
microphones currently available (including 
externally biased condenser, MEMS, and electret 
condenser devices), the authors chose to conduct 
a detailed performance analysis of a prototype 
DA system populated with COTS Electret 
Condenser Microphones (ECM’s). The 
objectives of the performance analysis were 
threefold: 

1. Characterize the acoustic response of a 
suitable ECM via laboratory testing and 
calibration, 

2. Evaluate the beamforming capability of an 
ECM-based directional array using a series 
of independently controlled point sources in 
an anechoic environment, and 

3. Demonstrate the utility of an ECM-based 
directional array in a real-world application. 
The authors chose as a demonstration to 
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apply the ECM array to measurement of 
source strength and directivity in a cold flow 
jet at high subsonic and low supersonic 
velocities. 

Note that the cost and performance of the 
sensors chosen to implement high-density DA 
systems are critical parameters determining the 
viability of such arrays, but they are not the only 
parameters to consider when designing such 
systems. Microphone package size, data 
handling efficiency (i.e., efficiency of the data 
acquisition system), and array mobility (i.e., 
number and size of cables, power connections, 
array mounting hardware, etc.) also determine 
the ultimate success of the system. While 
microphone package size and construction are 
discussed subsequently, array data handling and 
mobility are beyond the scope of this paper - 
reference 10 provides a more detailed discussion 
of these issues. Rather, it is the purpose of this 
paper to present an overview of the performance 
analysis conducted on the prototype ECM-based 
array constructed for this study. The authors 
believe this DA system represents a “small 
scale” version of what can be achieved using 
economical microphones. 

Microphone Construction 

The ECM’s employed for this work were 
“hybrid” devices composed of COTS 
microphone cartridges coupled with passive 
signal conditioning circuitry and packaged to 
enable easy integration into a DA system. The 
details of the packaging are described below. 

ECM Cartridge Several commercially-available 
ECM’s were investigated for use in construction 
of the prototype array system. After considering 
candidate devices, the authors chose to use 
Panasonic WM-60A ECM cartridges as the 
sensing element.* Figure la illustrates the 
construction of the ECM cartridge, which houses 
a recessed 0.25-inch diaphragm protected by a 
cap containing a central access hole 
approximately 0.13 inches in diameter. The 
microphone incorporates a field-effect transistor 


Specific manufacturer’s names are explicitly 
mentioned only to accurately describe the research 
performed. The use of manufacturer’s names does 
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nor does it imply that the specified equipment is the 
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driver on the output stage and can operate at a 
maximum supply voltage of +10V, with a 
maximum current draw of 0.5 mA. The nominal 
sensitivity of a raw WM-60A cartridge when 
driven by a +5V source through an 8.2 kQ load 
resistor is 15 mV/Pa. Figure lb illustrates the 
typical response of the microphones, which are 
designed to exhibit a flat magnitude response 
over the audio frequency range of 20 Hz - 20 
kHz. 

Equalizing Circuit / Packaging : While a typical 
WM-60A cartridge operates with a flat response 
over the audio frequency range, its use in 
aeroacoustic testing requires that the bandwidth 
of the device be extended to a significantly 
higher value. Above 20 kHz the microphone 
exhibits a first order sensitivity roll-off of 
approximately 6 dB per octave relative to its 
audio range sensitivity. To compensate for this 
roll-off, the authors designed a passive circuit to 
equalize the broadband frequency response of 
the microphone, in essence trading sensitivity for 
bandwidth. Figure 2 shows a schematic of the 
ECM cartridge with the equalizing circuit. The 
circuit was fabricated using surface mount 
devices on a small printed circuit board and 
mounted along with the ECM cartridge in a 
0.32-inch outside diameter aluminum cylinder 
approximately 2 inches in length. An 84-inch, 
four-conductor cable was attached to the rear of 
the cylinder via a gold-plated screw connector. 
The other end of the cable terminated in a quick 
disconnect coupler which was attached to a 
power and signal distribution patch panel. A 
total of 38 microphones were fabricated in this 
fashion, at an average cost per microphone of 
approximately $40. Figure 3 shows a 
representative finished microphone. 

ECM Calibration 

Four different types of calibration data were 
obtained - the 250-Hz sensitivity of each 
packaged microphone (raw cartridge coupled 
with equalizer), the broadband frequency 
response (magnitude and phase) of each 
packaged microphone, the background noise 
level for a sample of raw ECM cartridges, and 
the total harmonic distortion for a sample of raw 
cartridges. The various calibration procedures 
and representative data obtain from them are 
described subsequently. 

Sensitivity : Figure 4 shows a histogram of the 
250-Hz sensitivity distribution for the ensemble 
of 38 packaged microphones used in this study. 
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The sensitivity of each microphone was 
determined using a calibrated pistonphone 
generating a regulated Sound Pressure Level 
(SPL) of 94 dB (relative to a reference pressure 
of 20 pPa). The sensitivity of an individual 
microphone is defined here as the calibration 
value (specified in units of mV/Pa) which when 
applied to the output of the excited microphone 
allows it to indicate an SPL of exactly 94 dB at 
250 Hz. As can be seen in Figure 4, the 
sensitivities show a partial bimodal distribution 
centered around 6.25 and 7.75 mV/Pa, with a 
6.5 mV/Pa average sensitivity exhibited for the 
entire ensemble of microphones. The decrease 
in the sensitivity of the packaged microphones 
versus the raw WM-60A cartridges is due to the 
introduction of the equalizing circuitry. Note 
that the ensemble distribution is dominated by 
variations in the sensitivities of individual raw 
ECM cartridges only. Given the low-cost of 
these cartridges, if one desires a matched 
response among the microphones incorporated 
into a DA system, it is a simple matter to 
measure the sensitivities of a large sample of 
cartridges and then choose those with similar 
characteristics. However, this step was not 
performed for this study. Rather, the individual 
measured sensitivities of the packaged 
microphones were accounted for during analysis 
of all collected data. 

Broadband Frequency Response - Audio 
Spectrum : The broadband pressure response of 
each packaged microphone covering the audio 
frequency spectrum (250 Hz to 16000 Hz) was 
obtained using a standards-traceable Bruel and 
Kjaer 4226 multifunction calibrator which had 
been modified to allow access to the transducer 
driving signal. A special nylon coupler was 
fabricated to allow insertion of the microphones 
into the calibrator with minimal acoustic loss. 
The calibrator generated a constant SPL of 94 dB 
and was driven by an external precision signal 
generator operated over the range of 250 Hz - 
16 kHz in 250-Hz steps. The driving signal from 
the generator and the microphone output signal 
were recorded by a two-channel, PC-based 
transient data recorder with a per-channel 
sampling rate of 100 ksamples/sec. To obtain 
the magnitude and phase response of each 
microphone, the transfer function between the 
driving signal and microphone output were 
computed using Welch’s averaged periodogram 
method. 12 

Figure 5 shows the pressure responses 
obtained using the calibrator for 10 


representative microphones. The magnitude 
response of each device is referenced to its 
baseline 250-Hz sensitivity. As can be seen in 
the figure, the ensemble of magnitude responses 
(with one exception) matched to within ±5 dB at 
any one particular frequency while the phase 
responses matched to within ±10 degrees at a 
particular frequency over the range of 250 Hz 
through 16 kHz. The matching of phases within 
the range shown in Figure 5 implies that a 
directional array populated with these 

microphones should yield accurate beamformed 
estimates of noise strength without the need for a 
phase correction. The limits over which array 
phase matching needs to be performed are 
discussed at length by Mosher 13 . 

Broadband Frequency Response - Ultrasonic 
Spectrum : Magnitude and phase responses for 
each microphone for frequencies above the audio 
spectrum (16-40 kHz) were obtained using a 
substitution method freefield calibration. This 
calibration method conformed with the 

International Electrotechnical Commission (IEC) 
standard for calibration of microphone 
equipment. 14 Each fabricated electret 
microphone was mounted in an 8-foot by 12-foot 
by 6- foot high anechoic chamber directly in line 
with a series of 1.0-inch diameter ultrasonic 
ceramic transducers, as shown in the sketch in 
Figure 6. Each transducer was positioned 

74 inches from the microphone. The shortest 
microphone to chamber wedge-tip distance was 
approximately 24 inches. A custom mounting 
bracket with attached guide bar, shown in 
Figure 7, allowed correct alignment of each 
microphone in the chamber such that all 

diaphragms occupied as closely as possible the 
same axial location with respect to the source. 
Freefield calibration data were obtained at 
frequencies of 23.176, 29.976, and 39.576 kHz 
by operating each of three different ceramic 
transducers at its resonant frequency (to 
eliminate higher-order source diaphragm 
vibration modes). The transducer driving signal 
and the microphone output signal were recorded 
by the same two-channel PC-based transient data 
recorder used for the audio spectrum calibration, 
at a per-channel sampling rate of 
100 ksamples/sec. Upon completion of 
acquisition of calibration data using the three 
ceramic transducers, the electret microphone was 
replaced with a standards-traceable Bruel and 
Kjaer 4136 reference microphone and a repeat 
set of calibration data was acquired using the 
three transducers. The freefield magnitude and 
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phase response of the electret microphone at the 
resonant frequency of each transducer was 
obtained using Welch’s method in a similar 
manner to the audio spectrum calibration. 

One challenge encountered during the 
freefield calibration process concerned the 
determination of the proper corrections to apply 
to the freefield phase response. Such corrections 
were required to account for offsets introduced in 
the line of sight distance between the ceramic 
transducer and the microphone when installing 
different devices in the anechoic chamber, since 
the guide bar on the mounting bracket shown in 
Figure 7 provided a necessary, but not sufficient 
alignment. The preferred method for computing 
the phase correction was to generate a very short 
duration resonant frequency chirp using the 
ceramic transducer and then measure the time for 
the sound wave to travel from the source to each 
reference or test microphone. Unfortunately, the 
ceramic transducers employed for the test did not 
produce a chirp signature suitable for accurately 
determining the transit time. Thus, an alternate 
procedure was developed to obtain the phase 
corrections using a high frequency paper-cone 
tweeter in place of the ceramic transducers. The 
tweeter provided a much cleaner chirp signature 
for determining small transit time variations, and 
care was taken to not move the tweeter when 
installing different devices in the mounting 
fixture. Using an acquisition sampling rate of 
300 ksamples/sec, sufficient to resolve ~3 qsec 
time increments, the distance offset for each of 
the 38 microphones when installed into the 
anechoic chamber was adequately determined 
via a time series cross covariance analysis. Note 
that this process assumed that the reinstallation 
of a microphone into the mounting bracket using 
the guide bar resulted in an additional offset of 
no more than a few thousandths of an inch, 
corresponding to sub-microsecond travel time 
adjustments. As will be seen in the phase data 
shown subsequently, this correction procedure 
worked very well in practice. 

Figure 8 shows the freefield responses for ten 
representative packaged microphones obtained 
using the procedure described above. The 
freefield data (converted to pressure response) is 
overlaid with the Figure 5 audio spectrum 
pressure response data. As can be seen in the 
figure, the ultrasonic magnitude responses were 
matched to within ±5 dB at a given frequency 
while the phase responses were matched to 
within ±25 degrees. Since the high frequency 
phase response deviation between microphones 
greatly exceeded ±10 degrees, a phase correction 


was applied when using these devices in a DA 
system at frequencies above 20 kHz. Given that 
only three freefield calibration frequencies were 
acquired (due to the limited availability of 
ceramic transducers at ultrasonic frequencies 
below 40 kHz), a linear interpolation procedure 
was employed to provide the narrowband phase 
correction at frequencies other than the resonant 
ceramic ones. 

Background Noise Measurement : The 

background noise level of a raw ECM cartridge 
was determined using a vacuum-isolation 
chamber mounted on a vibration-damped table. 15 
The vacuum-isolation chamber consisted of an 
inner chamber maintained at ambient conditions 
into which the cartridge was placed. The inner 
chamber was then suspended within a vacuum 
chamber to provide a high degree of acoustic 
isolation from the surrounding environment. 
Figure 9 illustrates the background noise 
autospectrum for a typical cartridge for three 
different vacuum pressures. As can be seen in 
the figure, classic 1/f noise dominates the 
response at the lower end of the frequency 
spectrum. The background noise characteristic 
exhibited in Figure 9 is typical for electret-type 
microphones. 

Total Harmonic Distortion : The total harmonic 
distortion (THD) for a series of raw ECM 
cartridges was determined using an externally 
driven Whittaker PC- 125 acoustic calibrator 
excited at a frequency of 1 kHz for a range of 
sound pressure levels spanning 100 - 140 dB. 
The THD for a microphone is defined as the root 
mean square (RMS) value of the total harmonics 
of a measured calibration signal, divided by the 
RMS value of the fundamental plus harmonics of 
the signal. This can be expressed in percentage 
form as 

THD=^-x 100 % ( 1 ) 

P F 

where P H and P F are defined as 

P H = + P} ~\ b Z 3 ,, 2 (2a) 

P F = ^jp; + P{ + Pr +--- + P; (2b) 


5 

American Institute of Aeronautics and Astronautics 



measured over the first n harmonics. It was 
determined empirically that the fundamental and 
first three harmonics were sufficient to obtain an 
accurate estimate of the THD for the ECM 
cartridges tested. 

Figure 10 shows the THD measured on ten 
representative cartridges. Indicated on the plot is 
the average SPL level at which the THD reaches 
four percent. Sound pressure levels encountered 
during aeroacoustic testing in ground facilities 
(anechoic chambers, treated wind tunnels, etc.) 
typically do not exceed 120 dB; therefore, the 
maximum expected THD when using these 
microphones in a DA application is 
approximately 1-2%. 

Directional Array Evaluation System 

Array Construction : The performance analysis 
objective of testing the packaged electret 
microphones in a prototype DA was facilitated 
by utilizing an existing array mounting system. 
For this purpose, the authors chose to utilize the 
Langley Large Aperture Directional Array 
(LADA). 16 A photograph of the LADA is shown 
in Figure 11. A four-foot diameter fiberglass 
panel provided a flat surface to flush mount all 
microphones. The panel was attached to a pan- 
tilt unit secured to a rigid tripod support. This 
allowed precise alignment changes in the 
elevation and azimuth of the face of the array. A 
small laser diode pointer was place at x c , 
corresponding to the phase center of the array (in 
this case the center of the fiber glass panel), to 
aid in array alignment. The LADA incorporated 
35, 0.32-inch diameter microphone mounting 
holes arranged in a two-dimensional pattern 
consisting of logarithmic spirals. This 
microphone layout, shown in Figure 12, was 
based on a design by Underbrink et.al. at Boeing 
(see chapter 3 in reference 7). It consisted of 
five spirals of seven microphones each with the 
inner-most microphones lying on a 1-inch radius 
and the outer-most on a 17-inch radius. This 
design resulted in acceptable beamwidth and 
minimal sidelobe height over an operational 
frequency range of 2 - 30 kHz. The microphone 
mounting coordinates, referenced to x c , are 
listed in Table I. 

Array Population : The LADA was populated 

with either of two types of sensors - 
(1) commercial, standards-traceable, 0.25-inch 
externally-biased Bruel and Kjaer 4136 
condenser microphones (hereafter referred to as 


the “reference LADA”), and (2) fabricated and 
packaged electret microphones (hereafter 
referred to as the “ECM LADA”). The reference 
LADA microphones were powered by a bank of 
regulated DC supplies, one for each microphone. 
The ECM LADA was powered by a single 
regulated DC power supply, since the current 
draw from each electret microphone was no 
more than 0.3 mA using a supply voltage of 
+5 V and a load resistor of 8.2 kQ. In all cases 
the output signals from each microphone were 
conditioned using high-pass elliptic filters at a 
frequency of 300 Hz. Although not required, the 
electret microphone outputs were additionally 
amplified by a factor of 10 after filtering. 

Data Acquisition : The data acquisition system 
consisted of a 64-channel transient data recorder 
controlled by a PC. All data channels were 
simultaneously recorded with a 16-bit dynamic 
range at a sampling rate of 100 kHz. The analog 
bandwidth of the system at this sampling rate 
was 40 kHz based on the settings of the 
anti-aliasing filters. Typical acquisition run 
times ranged from 5 to 10 seconds with the data 
stored in raw binary format on the PC hard drive. 

Data Reduction : Post processing of acquired 

LADA data began with the computation of an 
M x M cross-spectral matrix for each ensemble 
of data, where M is the total number of 
microphones in the array. The computation of 
the individual matrix elements was performed 
using Fast Fourier Transforms (FFT) of the 
original data ensemble. This was done after 
converting the raw data to engineering units. The 
time data were segmented into a series of non- 
overlapping blocks each containing 4096 sam- 
ples. Using a Hamming window, each of these 
blocks of data was Fourier transformed into the 
frequency domain with a frequency resolution of 
24.4 Hz. The individual cross spectrum for 
microphones i and j was computed via 

G (/ (/) = ^7 X [4 (/)*,*(/)] (3a) 

where W s is the Hamming window weighting 
constant, N is the number of blocks of data, and 
X is an FFT data block. The full matrix was 
formed as 
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G = 


G u G n 

G 


22 


G 


1 M 


(3b) 


G 


MM . 


where the superscript T denotes a complex 
transpose of the matrix, and P(e ) is a mean- 
squared-pressure quantity. The division by M 2 
serves to reference the array output level to an 
equivalent single microphone output level. 
References (16) and (17) provide greater detail 
concerning additional processing performed on 
the LADA data. 


The lower triangular elements of this Hermitian 
matrix were computed by taking the complex 
conjugates of the upper triangular elements. All 
cross spectral matrix elements were employed in 
subsequent processing, with no modification of 
the diagonal terms. Note that for in-flow 
applications of the array, the diagonal terms can 
be removed to improve the spectral dynamic 
range by subtracting off self-noise dominated 
auto-spectra during the beamforming process 
(see chapter 1 of reference 7). However, this 
additional step was not required for the data 
acquired as part of this study. 

A classical delay and sum beamforming 
approach was used for the analysis of the array 
data. This approach assumed distributions of 
monopoles for the measured sources. In delay 
and sum beamforming, the array is electronically 
“steered” to a series of chosen source locations. 
For each selected steering location, a steering 
vector e containing a retarded time phase 
adjustment for each microphone in the array was 
defined as 


r 

-^expOXfc •*!>] 

r c 

— exp [j(k-xj] 

r 


where k is the acoustic wave vector, x m is the 
distance vector from the steering location to each 
microphone m, co is the angular frequency, and 
the ratio (r„/r c ) is included to normalize the 
distance related amplitude to that of the phase 
center of the array. The output power spectrum 
(or response) of the LADA at the steering 
location was obtained from 

/\T 

\ e Ue 

1 =^ 72 - ( 5 ) 



ECM LADA Evaluation 

Test Configuration : The beamforming capability 
of the ECM LADA was evaluated in the Langley 
Anechoic Noise Research Facility (ANRF). The 
ANRF is an anechoic chamber measuring 27.5 ft 
x 27.5 ft x 25.6 ft in which static or wind tunnel 
testing can be performed. The facility is 
primarily employed for engine nacelle noise 
studies. To assess the capability of the ECM 
LADA to identify discrete incoherent noise 
sources, independently controlled speakers were 
installed in the facility. These speakers were 
attached to mounting platforms positioned at 
various locations with respect to the center of the 
array. Figure 13 illustrates one such 
configuration where a single source (consisting 
of a piezo-electric cone tweeter) was located 
8 feet directly in front of the array. Both pure 
tone and white noise signals were recorded by 
the array microphones and processed using the 
data analysis procedures described previously. 
All calibrations performed on the ECM LADA 
were repeated for the reference LADA, using the 
same calibration frequencies, source strengths, 
and source locations. 

Performance Metrics : Two quantitative metrics 
were used to assess the accuracy of the 
beamformed output when using the ECM LADA 
to characterize point sources in the ANRF. The 
first metric consisted of computation of a 
weighted RMS error for the beamformed output 
of the ECM LADA compared with similar 
output from the reference LADA. Using an 
expression similar to that presented by Arnold, 
et.al. 18 , the weighted RMS error summed over all 
steering locations measured by the array can be 
expressed as 


£rms ~ 




fc = l 


K 


( 6 ) 
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where T k is the normalized measured response by 
the ECM LADA, R k is the normalized measured 
response by the reference LADA, and K is the 
total number of measured steering locations. 

The second performance metric consisted of 
comparing the beamformed output of the ECM 
LADA with its theoretical response. Assuming a 
simple monopole source, the LADA ideal array 
response can be expressed as 


W(k 9 x 9 x°) 



e M(r°-r)-(r°-r m )\ 


(7) 


where X is an arbitrary Cartesian location in 
space to which the array is electronically steered, 

x° is the source location, r° and r° are the 
distances from the source to X c and the m th 
microphone, respectively, and r and r m are the 
distances from the steering location to X c and the 
microphone, respectively. The response shown 
in eqn. (7) is normally expressed in decibels 

referenced to the level obtained at x° : 


dB(x ) = 201og 10 


W(k,x,x°) 

W(k,x°,x°) 


( 8 ) 


very similar between the theoretical response, 
ECM LADA output, and reference LADA 
output. The lower right portion of each figure 
shows a horizontal slice through each contour 
plot (normalized to the reference LADA output) 
and further illustrates the uniformity of the 
sidelobe structure when comparing the two 
arrays. 

Table II lists the weighted RMS errors for 
frequencies spanning the range of 
315 Hz -20 kHz, computed using eqn. (6) for 
normalized ECM- and reference-LADA 
beamform outputs similar to those shown in 
Figures 14 - 16. The error values are plotted in 
Figure 17. It is clearly seen in this figure that the 
weighted RMS errors are higher for very low 
beamform frequencies, but approach a 
frequency-invariant average value of 0.0046 
above 2500 Hz. This is a remarkable result and 
indicates that for the data obtained using a single 
point source in the ANRF, the beamform results 
when using the reference and ECM LADA’s are 
extremely consistent above a frequency of 
2500 Hz. Interestingly, this frequency also 
corresponds closely with the theoretical low 
frequency limit of the array microphone pattern 
for the LADA. 

Cold Flow Jet ECM LADA Demonstration 


Eqn. (8) is plotted as a color contour map with 
contour level proportional to dB(x) , computed 
over the complete ensemble of measured steering 
locations. Color contour maps generated using 
eqn. (8) over a range of beamform frequencies 
were compared with their experimental 

counterparts to assess the response fidelity of the 
ECM LADA. 

Sample Results : Figure 14 illustrates the 

theoretical and experimental array responses for 
the ECM and reference LADA’s measuring the 
5 -kHz narrowband component of a generated 
white noise source. The tweeter source (shown 
in Figure 13), was located with respect to the 
center of the LADA at (0.0, 0.0, 96.0) inches. 
The theoretical response was computed using 
eqns. (7) and (8), while the experimental 
beamformed output of the LADA was obtained 
using eqn. (5). Figures 15 and 16 illustrate 
similar responses for the 10-kHz and 20-kHz 
narrowband component of the white noise 
spectrum, respectively. As can be seen in these 
figures, the spatial distribution of sidelobes is 


Along with the beamforming evaluation using 
distributions of point sources, it was also desired 
to demonstrate the utility of the ECM LADA in a 
real-world application. The authors chose as a 
demonstration to apply the array in a series of 
simple measurements of source strength and 
directivity in a cold flow jet operated at high 
subsonic exit velocities in the ANRF. Note that 
it is beyond the scope of this paper to present a 
detailed study of jet noise using the array. 
Rather, the goal of the demonstration was to 
apply and evaluate the ECM LADA in a rich 
sound field containing broadband noise and high 
background levels. 

The jet which was chosen for the 
demonstration was a 2-inch diameter, Mach 1.0 
contour convergent nozzle which had been 
previously used for measurements of mean 
pressure and farfield acoustics. 19 Figure 18 
shows a photo of the nozzle mounted in the 
ANRF with the ECM LADA shown in the 
background. The majority of the nozzle shown in 
the photo was covered with acoustic foam to 
prevent sound reflections. The jet was operated 
at high subsonic exit velocities spanning Mach 
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numbers from 0.50 to 0.99 to provide generated 
noise which was then measured by the array. 
Noise data were acquired at several orientations 
of the array with respect to the jet exit plane (i.e., 
various azimuth angles with respect to the jet 
axis) to obtain a measure of the jet directivity. 
Figure 19 shows a sample noise location map for 
a narrowband frequency of 5 kHz obtained with 
the LADA positioned 8 feet from the center of 
the nozzle, pointing directly at the exit plane 
(azimuth angle = 90 degrees). 

One useful characteristic of a DA system 
when measuring free jet flow noise is the array’s 
ability to extract from the beamformed data the 
variation of the axial location of the jet peak 
source strength as a function of frequency. As 
discussed by Lilley 20 , the axial location data 
should collapse to single curve when the 
frequency parameter is expressed as the Strouhal 
number 


Panasonic WM-60A electret cartridge were 
fabricated for use with the Langley Large 
Aperture Directional Array. The results of the 
investigation revealed a fabricated microphone 
frequency response suitable for directional array 
use, with minimal magnitude and phase response 
variations in the audio frequency spectrum and 
moderate variations at higher frequencies. Also, 
the results using an ECM-populated LADA to 
measure point source noise in the Langley 
ANRF revealed excellent agreement between 
reference array data, theoretical array response, 
and the LADA beamformed output. Finally, a 
successful demonstration was conducted using 
the ECM LADA to measure properties of the 
noise field generated by a cold flow jet. Thus, 
economical COTS electret microphones show 
great promise for array work, especially for 
lower-frequency applications such as the jet 
noise demonstration presented here. Additional 
studies similar to this one are planned for COTS 
MEMS microphones which have recently 
become available. 


where / is the acoustic frequency, Dj is the jet 
diameter, and Uj is the jet exit velocity. 
Figure 20 shows the axial location of the jet peak 
source strength as measured by the ECM LADA 
for three different jet exit velocities spanning 
Mach numbers from 0.5 to 0.99. Also shown in 
the figure are two theoretical curves presented by 
Lilley in reference 20 which were obtained from 
calculations based on Lighthill’s acoustic 
analogy. The lower curve is based on the 

assumption of a zero initial thickness for the jet 
mixing region. The upper curve is based on the 
assumption that the mixing region becomes fully 
turbulent a defined distance downstream from 
the jet exit. These two curves bracket the ECM 
LADA results, indicating a good ability of the 
array to extract useful information from the jet 
noise data. 


Summary 


Directional array technology will continue to 
be an integral part of experimental aeroacoustics. 
However, cost remains a limiting factor in the 
construction of future high-density array systems 
which will be needed to implement advanced 
beamforming algorithms. In an effort to reduce 
the construction costs of future arrays, the 
authors have undertaken a systematic 
investigation of one class of economical sensing 
device, the electret condenser microphone. A 
series of customized microphones based on the 
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Table I - LADA Microphone Coordinates (Viewed from Front of Array) 


Mic# 

X 

location 

Y 

location 

Z 

location 

Mic# 

X 

location 

Y 

location 

Z 

location 


(in) 

(in) 

(in) 


(in) 

(in) 

(in) 

1 

0.03 

-1.02 

0.00 

19 

-9.99 

5.66 

0.00 

2 

0.97 

-0.32 

0.00 

20 

-8.47 

-7.78 

0.00 

3 

0.61 

0.80 

0.00 

21 

10.57 

-8.60 

0.00 

4 

-0.57 

0.81 

0.00 

22 

11.44 

7.37 

0.00 

5 

-0.93 

-0.30 

0.00 

23 

-3.47 

13.14 

0.00 

6 

-5.13 

0.88 

0.00 

24 

-13.56 

0.75 

0.00 

7 

-2.42 

-4.63 

0.00 

25 

-4.89 

-12.69 

0.00 

8 

3.67 

-3.75 

0.00 

26 

14.55 

-5.15 

0.00 

9 

4.71 

2.30 

0.00 

27 

9.39 

12.22 

0.00 

10 

-0.73 

5.17 

0.00 

28 

-8.73 

12.68 

0.00 

11 

-2.34 

-8.63 

0.00 

29 

-14.74 

-4.38 

0.00 

12 

7.50 

-4.91 

0.00 

30 

-0.39 

-15.42 

0.00 

13 

6.99 

5.59 

0.00 

31 

17.01 

-1.12 

0.00 

14 

-3.16 

8.35 

0.00 

32 

6.32 

15.81 

0.00 

15 

-8.92 

-0.42 

0.00 

33 

-13.08 

10.89 

0.00 

16 

4.80 

-10.49 

0.00 

34 

-14.39 

-9.09 

0.00 

17 

11.45 

1.30 

0.00 

35 

4.21 

-16.52 

0.00 

18 

2.31 

11.28 

0.00 






Table II - RMS Errors for ECM LADA with Tonal Point Source 


Frequency, Hz 

Srms 

315.0 

0.0591 

400.0 

0.0288 

500.0 

0.0409 

630.0 

0.0294 

800.0 

0.0279 

1000.0 

0.0112 

1250.0 

0.0134 

1600.0 

0.0146 

2000.0 

0.0102 

2500.0 

0.0047 

3150.0 

0.0041 

4000.0 

0.0060 

5000.0 

0.0051 

6300.0 

0.0040 

8000.0 

0.0036 

10000.0 

0.0037 

12500.0 

0.0038 

16000.0 

0.0056 

20000.0 

0.0058 
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Figure 1. (a) Internal Construction of Electret Microphone, 
(b) Manufacturer’s Audio Frequency Response. 
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Figure 2. Schematic of Electret Microphone, Power Loading, and Passive Equalizer. 



Figure 3. Representative Fabricated Microphone and Cable. 
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Phase, deg dB re 250 Hz Sensitivity 



Sensitivity, mV/Pa 


Figure 4. Histogram of Sensitivities for Packaged Microphones. 




Figure 5. Audio Spectrum Response of 10 Representative Packaged Microphones. 
Measured Using B&K 4226 Multifunction Calibrator. 
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Not to Scale 


Figure 6. Sketch of Anechoic Chamber Figure 7. Photo of Freefield 

Freefield Calibration Configuration. Mounting Bracket with Attached Guide Bar 
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Figure 8. Total Response of 10 Representative Packaged Microphones. 
Left Side - Audio Spectrum Calibration using B&K 4226 Calibrator. 
Right Side - Freefield Ultrasonic Transducer Calibration. 
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Figure 9. Typical Vacuum-Isolation Chamber Background Spectra 
for Raw ECM Cartridge. 
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Figure 10. Total Harmonic Distortion of Ensemble of Raw ECM Cartridges. 
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Figure 11. Langley Large Aperture Directional Array. 
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Figure 13. ANRF Point Source / LADA Configuration. Point Source is 
Located at Lower Left of Image, LADA at Center Right. 
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Figure 14. LADA Point Source Run in LaRC ANRF. Source at (0.0, 0.0, 96.0)” wrt Array Center. 
Upper Left - Theoretical Sidelobe Response for 5-kHz Signal. 

Upper Right - Reference LADA Response. 

Lower Left - ECM LADA Response. 

Lower Right - Horizontal Slice Through Responses (red=theory, green=reference, blue=ECM). 
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Figure 15. LADA Point Source Run in LaRC ANRF. Source at (0.0, 0.0, 96.0)” wrt Array Center. 
Upper Left - Theoretical Sidelobe Response for 10-kHz Signal. 

Upper Right - Reference LADA Response. 

Lower Left - ECM LADA Response. 

Lower Right - Horizontal Slice Through Responses (red=theory, green=reference, blue=ECM). 
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Figure 16. LADA Point Source Run in LaRC ANRF. Source at (0.0, 0.0, 96.0)” wrt Array Center. 
Upper Left - Theoretical Sidelobe Response for 20-kHz Signal. 

Upper Right - Reference LADA Response. 

Lower Left - ECM LADA Response. 

Lower Right - Horizontal Slice Through Responses (red=theory, green=reference, blue=ECM). 
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Figure 17. Weighted-RMS Errors Between ECM and Reference LADA 
Beamformed Output for Single Point Source in LaRC ANRF. 



Figure 18. Acoustically Treated Nozzle and ECM LADA in ANRF. 
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Figure 19. Typical Jet Flow ECM LADA Source Map. M=0.647, Frequency=5000 Hz. 
Steering Grid Overlaid on Map. Dj=2 inches. 
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Figure 20. Axial Location of Jet Peak Source Strength for Range of Exit Velocities. 
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